Quantum Algorithm for Solving an NP-Complete Problem 
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When a probe qubit is coupled to a quantum register that represents a physical system, the 
probe qubit will exhibit a dynamical response only when it is resonant with a transition in the 
system. Using this principle, we propose a quantum algorithm for solving a specific NP-complete 
problem, the 3-bit Exact Cover problem, EC3. We show that on a quantum computer, the number 
of qubits increases linearly with the size of the EC3 problem, while the efficiency of the algorithm is 
independent of the size of the problem. Our results indicate that quantum computers may be able 
to outperform classical computers in solving NP-complete problems. 

PACS numbers: 03.67.Ac, 03.67.Lx 



Extraordinary progress has been made in the field of 
quantum computation since the discovery of Shor's fac- 
toring algorithm [l[ and Grover's search algorithm Q. It 
was found that quantum computing offers speedup for 
a number of problems 0, Classically, problems can 
be divided into two groups, those that can be solved in 
polynomial time (P) and the nondeterministic polyno- 
mial (NP) problems for which no polynomial-time algo- 
rithms are known. NP-complete (NPC) problems are a 
subclass of NP problems. Solving one NPC problem in 
polynomial time means that all NPC problems can be 
solved in polynomial time. The question of whether there 
are polynomial-time algorithms for NPC problems is an 
unsolved problem in mathematics [f|. It is also not 
known whether or not NPC problems can be solved in 
polynomial time on quantum computers. In this paper, 
using the principle that when a probe qubit is coupled to 
a quantum register that represents a physical system the 
probe qubit will exhibit a dynamical response only when 
it is resonant with a transition in the system, we pro- 
pose a polynomial-time quantum algorithm for solving a 
specific NPC problem, the 3-bit Exact Cover problem, 
EC3. 

The EC3 problem is a restricted version of the 3-bit 
satisfiability (3-SAT) problem Q. No algorithm has been 
found that can solve this problem in polynomial time. In 
Ref. @, Hj], Farhi and coworkers developed a quantum 
adiabatic algorithm (QAA) for solving the EC3 problem. 
In the QAA for solving the EC3 problem, one starts from 
an initial Hamiltonian, Hs, and its ground state as the 
initial state, and evolves Hb to a final Hamiltonian Hp, 
whose ground state is the solution to the EC3 problem. 
The system is evolved from the initial state to the ground 
state of Hp. The computation time t increases as t ~ 
Smin> where <?min is the minimum energy gap between the 
ground and the first excited state during the evolution. 
The cost of the QAA is determined by the scaling of 
Smin with the size of the EC3 problem (here, the size of 
the EC3 problem is defined as the number of bits in the 
problem). It has been found that the QAA cannot solve 
the 3-SAT problem, or the EC3 problem, efficiently [§- 
13j . Here by applying a quantum algorithm we developed 



recently |14| , we present a quantum algorithm for solving 
the EC3 problem. 

The EC3 problem on a quantum computer can be for- 
mulated as follows 0, Q : the 3-bit instance of satisfiabil- 
ity is a Boolean formula with M clauses 
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where each clause Ci is true or false depending on the 
values of some subset of the n bits, and each clause con- 
tains three bits. The clause is true if and only if one of 
the three bits is 1 and the other two are 0. Let ic, jc 
and kc be the 3 bits associated with clause C. For each 
clause C define an "energy" function 



Then 
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Jo, if (zi c , Zj c , Zk c ) satisfies clause C 
|l, if (zi c , zj c , Zk c ) violates clause C 

(2) 



Hp,c\z\Z2 ■■■Z n ) =hc (Zi ,Zj G ,Zk ) \z\Zi ■ ■ ■ z n ), (3) 
where \zj) is the j-th bit and has the value or 1. Define 



(4) 



Hp\if>) = 0, if and only if \tjj} is a superposition of states 
of the form \z1Z2 ■ ■ ■ z n ), where the bit string Z\Z<i ■ ■ ■ z n 
satisfies all of the clauses. If formula (1) has no satis- 
fying assignments, the ground state (or states) of Hp 
corresponds to the assignment (or assignments) that vio- 
lates the smallest number of clauses. The computational 
basis of \z1Z2 ■ ■ • z n ) is of dimension N = 2™. The task is 
to find out if there is one (or more) of the 2 n assignments 
that satisfies all of the clauses, that is, makes formula (1) 
true. 

In the previous work we proposed a quantum al- 
gorithm for obtaining the energy spectrum of a physical 
system. In this algorithm, we couple the system to a 
probe qubit with a certain frequency, set the probe qubit 
in its excited state, evolve the whole system for some 
time, then perform a measurement on the probe qubit. 
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When the frequency of the probe qubit matches the tran- 
sition frequency between two energy levels of the system, 
one observes a peak in the decay rate of the probe qubit. 
Therefore by varying the frequency of the probe qubit, 
we can locate the transition frequencies of the system. 
With a little variation, this algorithm can also be used 
for solving the EC3 problem. The details are shown be- 
low. 

We construct a Hamiltonian H that has the form: 



H 



-In 
Hp 



(5) 



where In is the iV-dimensional identity operator, and H 
acts on the state space of a (n+l)-qubit quantum register 
Rs- We let a probe qubit couple to Rs, and design a 
Hamiltonian H for the whole system of the form 
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where the first term is the Hamiltonian of the probe 
qubit, the second term is the Hamiltonian for the quan- 
tum register Rs, and the third term describes the inter- 
action between the probe qubit and Rs- Here, u is the 
frequency of the probe qubit (we have set K = 1), and 
c is the coupling strength between the probe qubit and 
Rs, while a x and a z are Pauli matrices. The operator 
A acts on the state space of H and plays the role of an 
excitation operator. It contains all possible excitation 
operators that excites the states from subspace of —In 
to subspace of Hp, 



A=-L(A 1 + A 2 + 
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(7) 
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where Ai = a x ® I 2 ® h <£> • • • ® h, A 2 = a x 
■ ■ ■ <g> I 2 (g> <7 X , • • • , An+i = a x (8 a x ® I 2 65 I 2 ® • 
A n+2 — a x ® 1 2 ®- ■ - ® 1 2 ® a X ®(J X , A n+ 3 = tr x ®/a(8- • -® 
<j x ®I 2 ®cf x , . . . , An = (J x (g)(T x (g)- ■ ■<g)(T x (g)(T x , where I 2 is 
the two-dimensional identity operator. Even though the 
operator A has an exponentially large number of terms, 
it can be implemented at polynomial cost as 
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We prepare Rs in the state 



JV N 



(8) 



(9) 



as the reference state, where IxJ are the states of the 
computational basis in the subspace of In, and |0) = 

q ^ . This is achieved by initializing Rs in state 

|0)®("+i) and applying the operator I 2 <E> Hf n , where 
Hd is the Hadamard gate. The states |<Pj) are the eigen- 
states of H with eigenvalues —1 and with a degeneracy of 
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FIG. 1: Quantum circuit for solving the Exact Cover problem 
EC3. The first line represents a probe qubit, and the next n+1 
lines represents the quantum register Rs- Hd represents the 
Hadamard gate, and U(t) is the time evolution operator with 
Hamiltonian given in Eq. (6). 



N, therefore the reference state |\& s ) has the eigenvalue 
— 1. We set the frequency of the probe qubit u> = 1, and 
run the circuit in Fig. 1. Let the whole system evolve 
with the Hamiltonian shown in Eq. (6) for time r. After 
that, we read out the state of the probe qubit. We repeat 
this procedure many times in order to obtain the decay 
probability of the probe qubit. 

In the procedure above, setting the probe qubit fre- 
quency u> = 1, we can probe the transition between the 
initial state, whose eigenvalue is — 1 , and the eigenstate 
of H with corresponding eigenvalue 0, which encodes the 
solution to the EC3 problem. If there exists at least one 
assignment that satisfies all of the clauses, a decay of the 
probe qubit will be observed. The assignment is encoded 
in the state of the last n qubits of the quantum regis- 
ter Rs- If there is no assignment that satisfies all of the 
clauses, we can increase a; by 1 each time and run the 
algorithm. Since the eigenvalues of the Hamiltonian Hp 
are integers, if we keep increasing the frequency of the 
probe qubit by 1 and run the algorithm, we will find the 
assignment that violates the smallest number of clauses. 

The procedure of the algorithm is summarized as fol- 
lows: (i) prepare a probe qubit in state |1) and a (n+1)- 
qubit quantum register Rs in state |0)®(™ +1 ); (ii) apply 
operator I 2 <8> H? n on Rs, therefore Rs is transformed 
into state \^ s ) as shown in Eq. (9); (Hi) implement the 
time evolution operator U(t) = cxp (— iifr) where H is 
given in Eq. (6); (iv) read out the state of the probe 
qubit; (v) perform steps (i) - (iv) many times in order to 
obtain good statistics and calculate the decay probability 
of the probe qubit. The quantum circuit for steps (i) — 
(iv) is shown in Fig. 1. 

In our algorithm, n + 2 qubits are required for solving 
an n-bit EC3 problem on a quantum computer, which 
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increases linearly with the size of the problem. In the 
following, we discuss the efficiency of the algorithm. 

The efficiency of the algorithm is defined as the number 
of times that the circuit in Fig. 1 must be run in order 
to obtain the decay probability, Pdecay The number of 
times that the circuit needs to be run must be at least 
proportional to 1/Pjccay in order to obtain Pdecay 

As we discussed in Ref. the decay probability of 
the probe qubit is: 



decay " 



where Qi 



2 f*^ij T 



Qf+iEj-Ei-J) 



d(A,|* s )| 2 , (10) 



2c|(A j |A|A i )|, and 



we have weak system-probe coupling and low noise. The 
evolution time r should be large, such that the change 
of the system is clear and one obtains high decay prob- 
ability. From Eq. (11), we can see that the efficiency 
of the algorithm is independent of the size of the prob- 
lem. Therefore by choosing appropriate c and r, the EC3 
problem can be solved efficiently on a quantum computer 
using our algorithm. 

When we set uj = 1, the coupling between the initial 
state and all the other states except the state with eigen- 
value equals to zero also contribute to the decay proba- 
bility of the probe qubit, therefore introduce an error in 
■Pdecay We now evaluate this error, Pdecay 
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Qij + {E 3 - E t - u) . \Xi) is the i-th energy pclI _ V^ sin 
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eigenstate and E t is the corresponding eigenenergy of 
the system, and according to Eq. (4), P, ; are discrete 
integers. Eq. (10) describes Rabi-oscillation dynamics, 
where the system and probe exchange an excitation. 

If a solution to the EC3 problem exists, the excitation 
frequency between the initial state |\E' S ) and the state 
with eigenvalue 0, which is the solution to the problem, 
is 1. With the probe qubit frequency being set to u> = 1, 
and considering the degeneracy of !<£>,), the decay proba- 
bility of the probe qubit becomes 



P 



decay 



= Sill 



2 



(11) 



where 



O 



2cMA\* s )\ 

N m , 

i=i j=i v 

2cy/m. (12) 



This term describes the summation over all excita- 
tion channels in the Rabi-oscillation. Here = 
jyjLi |l)lMj)/'V / ^ encodes the solution to the EC3 prob- 
lem which is a superposition of all m assignments that 
satisfy all of the clauses. And |// •) are these assignments, 
which are the basis states of Hp with eigenvalues and 
have degeneracy m. From Eq. (12), we can see that as 
m increases, there will be more excitation channels and 
the period of Pdecay will decrease. By knowing fio, we 
can choose a specific evolution time r by considering the 
degeneracy of the ground state of the problem, such that 
the probe qubit has high decay probability, therefore the 
algorithm has high efficiency. By trying a few times, we 
can find the optimal evolution time r for implementing 
the algorithm. 

From the equations above, we can see that the decay 
probability, and therefore the efficiency of the algorithm, 
depends on the coupling strength c and the evolution 
time r. In general, we need to set c to be small so that 
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where Eq — —1, is the eigenvalue of the reference state, 
TTtj represents the degeneracy of the assignments that 
have eigenvalue Ej. N m is the eigenenergy of the high- 
est energy level of the problem, and as N m — > oo, 
mo is the maximum of m,-. The 



term 2mo7r 2 c 2 /3 can be small since c can be set small. 
Therefore we can control this error and make it very 
small. 

In the following, we present an example that demon- 
strates how the algorithm would perform. 

Consider an 8-bit EC3 problem, we choose three cases 
of the 3-bit sets of clauses, which have one, two and four 
satisfying assignments to the problem, respectively. As 
discussed above, 10 qubits are required for solving this 
problem using our algorithm. We set the coupling coef- 
ficient c = 0.002, and run the algorithm with different 
evolution time r. As the probe qubit decays to the state 
|0), the state of the last 8 qubits of the quantum register 
Rs encodes the solution to the above EC3 problem. Case 
i): the 3-bit sets are {1, 2, 8}, {2, 3, 6}, {2, 3, 7}, {2, 4, 5}, 
{2,5,6}, and {3,5,8}. The solution to the 8-bit EC3 
problem for this case is |00010111). Case ii): the 3-bit 
sets are chosen to be {1, 4, 5}, {1, 7, 8}, {2, 4, 8}, {2, 7, 8}, 
{4,5,8}, and {5,6,7}. The solution to the problem for 
this case is |00010010) and |00110010). Case Hi): the 
3-bit sets are {1,3,5}, {1,6,8}, {2,4,6}, {2,6,8} and 
{4,5,7}. The solution to the problem for this case is 
100001100), 100100110), 100110001), and |11000010). 

In Fig. 2, we show the variation of simulated Pdecay 
with the evolution time r for the above three cases. The 



4 




400 800 1200 1600 

Evolution time r 



FIG. 2: Decay probability of the probe qubit vs. evolution 
time t. The frequency of the probe qubit u; = 1 and the 
coupling strength c = 0.002. The black solid line shows the 
result for case i): the 3-bit sets are {1, 2, 8}, {2, 3, 6}, {2, 3, 7}, 
{2,4,5}, {2,5,6}, and {3,5,8}; the red dashed line shows 
the result for case ii): the 3-bit sets are {1,4,5}, {1,7,8}, 
{2, 4, 8}, {2, 7, 8}, {4, 5, 8}, and {5, 6, 7}; and the blue dotted 
line shows the result for case Hi): the 3-bit sets are {1, 3, 5}, 
{1,6,8}, {2,4,6}, {2,6,8} and {4,5,7}. 



black solid line, the red dashed line and the blue dotted 
line show the results for cases i), ii) and Hi), respec- 
tively. In the above three cases, the decay probability 
almost reaches one at r = 800, r = 550 and r = 400, 
respectively. This shows that the algorithm can be run 
efficiently From the figure we can also see that as the 
number of satisfying assignments increases, the period of 
-Pdecay decreases. The simulated Pdecay fits well to the 
analytical results predicted by Eq. (11). 

We now discuss the implementation of the algorithm. 
In the algorithm, we have to implement the time evolu- 
tion operator U(t) = exp (— iHr). In the Hamiltonian 
H, the first two terms commute, while they do not com- 
mute with the third term. As shown in Eq. (8), the third 
term can be implemented efficiently using only (n + 2) 
gates. The operator U(t) can be implemented using the 
procedure of quantum simulation based on the Trotter- 
Suzuki formula 11511: 



U(r) 



\o(\), (14) 



where L does not depend on the size and the Hamiltonian 
of the problem, it can always be chosen sufficiently large 
such that the error is bounded by some threshold [la ]. 
The second term of the Hamiltonian H , H, as shown in 
Eq. (5), can be seen as a Controlled- Hp (C-Hp) opera- 
tion. The operation that calculates Hp can be taken as 
an oracle. In our algorithm, this oracle is entangled with 
the probe qubit. The number of times that the C-Hp or- 
acle is implemented is L. Therefore the implementation 
of U(t) scales polynomially with the size of the problem. 
In our algorithm for solving the EC3 problem, we con- 



struct a Hamiltonian that contains the Hamiltonian of 
the problem and a Hamiltonian with the same dimension 
N of the problem, whose eigenstates have degeneracy N 
and eigenvalues (in our case, —1) lower than the small- 
est eigenvalue of the problem, and are used as a refer- 
ence point. If a solution to the problem exists, we will 
observe the decay of the probe qubit at the excitation 
frequency (in our case, 1) between the reference state 
and the eigenstates with eigenvalue zero. In this case the 
system register evolves to its ground state, which is a su- 
perposition state of the assignments to the problem. If 
there is no solution to the problem, we can increase the 
frequency of the probe qubit discretely (since the energy 
function of the problem is discrete) , and record the decay 
probabilities of the probe qubit. The first frequency at 
which the decay of the probe qubit indicates the ground 
state of the system, which is the solution that violates the 
smallest number of clauses. This algorithm can also be 
used for efficiently solving some other discrete mathemat- 
ical problems, such as some combinatorial optimization 
problems. 
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